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Abstract 

An extension of the Zee model involving a light right handed neutrino, vr is con- 
sidered. We update constraints on couplings between the bilepton scalar, the active 
neutrinos, vr and the charged leptons. We find that the most stringent constraint 
currently comes from measurements limiting the width of the decay fj, — ► e 7. These 
are used to predict the upper bound on violation of lepton universality in leptonic 
W boson decays and rare Z decays, such as Z — ► e/x. 

PACS: 13.15+g, 12.60-i 



Evidence is now mounting that points towards neutrinos having small but finite 
masses. The neutrino mixing explanations which have been invoked to explain the deficit 
of solar neutrinos, the atmospheric neutrino anomaly and the LSND results, require neu- 
trino mass squared differences in the range, Sm 2 of 10~ 6 to a few eV 2 . The exact values 
depend on the details of the models used to analyze the data. The standard explanations 
for each phenomenon involve mixing between two flavors of neutrino. For example the 
atmospheric neutrino anomaly is considered to be large angle mixing between the muon 
neutrino and something else, which is not the electron neutrino. These three standard 
interpretations are not easily reconciled without four neutrino states. In an alternative, 
three neutrino mixing, it is difficult to accommodate the zenith angle dependence of the 
Super Kamiokande atmospheric results 0. 

It is generally accepted that neutrinos with masses in the above mentioned range 
necessitate the extension of the Standard model (SM) of electroweak interactions. Early 
simple models involving neutrino mass and/or a fourth neutrino state were constructed in 
and J|. The simplest way to generate neutrino mass is to add at least one SM singlet 
right- handed fermion, denoted by Ur, to the matter content. Indeed Grand unified models 
such as SO (10) naturally contain such a neutral fermion. In such a model, masses for the 
ordinary neutrinos can be elegantly explained by the seesaw mechanism. For masses as 
small as those predicted by solar, atmospheric and LSND neutrino mixing, the SM singlet, 
vr is required to have a mass in the range of 10 10 — 10 12 GeV. The seesaw mechanism does 
not easily accommodate a light fourth neutrino. Therefore, in accepting such a scheme 
one also implicitly accepts the view that one or more of the explanations of the evidence 
for neutrino mixing is misleading. Furthermore, there are several uncertainties which arise 
when utilizing the seesaw mechanism to produce light neutrinos. The means of obtaining 
this high energy SM singlet is highly model dependent. The energy scale of the heavy vr 



1 



is arbitrary and reproducing any of the mass and mixing schemes usually requires model 
embellishments. 

In view of these uncertainties it is important to investigate scenarios or models which 
relate the small neutrino mass with new physics at the weak scale or just above, i.e. 
< O(TeV). This alternative has the phenomeno logical advantage that such models may 
be easily constrained or perhaps tested. The simplest scenario was constructed by Zee|J 
many years ago. In the original formulation only the ordinary left-handed neutrinos 
are employed and their masses are generated by radiative corrections. Since the masses 
produced in this way are naturally small, it is unnecessary to invoke a large mass scale. A 
light vr can also be incorporated as an extension to the model. However, this is done at 
the price of giving up the predictability of the neutrino masses and mixings || [7|, || . To 
date the bulk of the literature on the Zee model is devoted to the study of the neutrino 
mass matrix and neutrino mixings. The charged scalar meson, S~ , which we shall refer 
to as the bilepton scalar is the crucial agent that carries the lepton flavor violation (LFV) 
necessary for neutrino mass generation and has been relegated to a secondary role in all 
studies thus far. 

In this paper we study the Zee model as the simplest model of lepton number violation 
at the weak scale. In view of all of the data on neutrino oscillations, we have augmented it 
with a vr [0] . We focus on the physics involved with virtual effects of the charged scalar. 
Since no charged scalar has been found up to LEPII energies, we assume its mass to be 
greater than 100 GeV. Due to the simplicity of this model, the number of free parameters 
is relatively few. Furthermore, these can be tightly constrained by current precision 
measurements. We begin by updating all the constraints stemming from muon decays 
and tau decays. We also discuss the impact of the next generation of Michel parameters 
measurements |10| on the model. We find that the strongest constraint currently comes 
from jj, — > decay. We then present new results for leptonic universality tests involving 
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physical W boson decays. This is particularly important in view of the large sample of 
W bosons which will be obtainable from the LHC and the next linear collider (NLC). 
We also calculate corrections to the left-handed charged lepton and Z boson couplings as 
well as the right-handed charged lepton and Z couplings. As we shall see below these two 
types of corrections take very different forms. We also give the predicted widths of the 
rare Z decays such as Z — > e/x. As far as we can determine these results have not been 
presented before. 

Without further ado the following Lagrangian is added to the SM: 

£ s =[/i2&L - ^e L ) + f 13 (v°T L - T*e L ) + / 23 (^t l - v*n L ) 

_ _ _ C 1 ) 

+ g\v%e R + gzv%\i R + g 3 u c R r R ]S + + h.c. , 

where S + is the Zee scalar boson with hypercharge Y=2. The and gi are Yukawa 
couplings and make up six free coupling parameters of the model. If the all the g's are 
set to zero one recovers the original Zee model. The SU(2)xU(l) charges and the lepton 
number L of the leptons and S~ are presented in Table I. It is seen that the Lagrangian 
conserves total lepton number; however, individual e, mu or tau number is violated. 
We have not included the scalar potential for the S~ and its interaction with the SM 
Higgs doublet as they are not needed here. It suffices to note that the Zee boson can 
develop a mass either by spontaneous symmetry breaking via coupling to the SM Higgs 
doublet and/or explicitly through a bare mass term. In either case the physical mass is 
another free parameter which we denote by M$- The above Lagrangian has been used to 
study neutrino oscillations in ||. (It was found that the v R significantly changes the 
phenomenology of neutrino oscillations from that of the three light neutrinos case and a 
fit to all data can be achieved). It will be seen below that since we have assumed a light 
mass for this particle, it can impact precision electroweak measurements in unusual ways, 
due to its chirality. Significant bounds on the interaction strength of the v R can already 
be obtained with presently available data. 
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Since LEP II has set Ms to be higher than the W mass at low energies we can integrate 
out the S boson and perform a Fierz transformation to obtain the following effective four- 
fermion Lagrangian in terms of weak eigenstates: 

-£l ff = 2M^\ l^ 12 ' 2 ^ 7 ' 3 ^ 7 ^ + f 12 f 13^11 p^Y^l (2) 
+9*lf eL^R^L ~ 92 f^J^^P-R ~ 92fv^lVRl\PR 

+1/ l 4:[glf2 3 e R -o p \v T v^a p> " p L + g 2 flz^L° ' P \ V R^ P> 'PR 

+92 'p\^R~p(y p) ^PR ~ g\f\2^<5p\VeVRO pX PL\^ + tl.C. . 

We have shown only the terms relevant for \x decay. The dominant term for the muon 
lifetime is then given by 

AG SM 

~ £eff = —J=~ {gZljLlpVeLVpJ.YPL + •••}, (3) 

where 

gl h = 1 + Ag LL 

2\f 12 \ 2 M* w (4) 
9 2 Ml ' 

and g is the SU(2) gauge coupling. 

We have adhered to the notations of the Particle Data Group jTT) . For the purpose of 
this paper it is convenient to work in the weak eigenbasis. One can easily see from Eq(0) 
that without z/r the Zee model will give rise the same chiral structure as the SM [ see 
Eq.(BI)]. The first term will interfere coherently with the SM amplitude and is the most 
important contribution to the muon lifetime, r M . The terms involving Vr add incoherently 



can be neglected in the muon lifetime r^. Explicitly, [JT2]1|]I3" | 



where we have used the radiatively corrected expression and have neglected terms involv- 
ing m e /M w and m^/M^. The SM fermi constant is given by ||14j| 



qSM = ™ /g\ 

F V2M W (1 - M^/Mf) (1 - Ar) ' 1 ' 

Since is one of the most accurately determined quantity in particle physics a careful 
analysis is need in order to obtain a bound on the coupling versus Ms. In Eq. (|6]) 
the fine structure constant a(0) is very accurately known to be 1/137.03598959 and Ar 
is the SM correction and is found to be ]T3| Ar = .0349 =F .0019 ± .0007. However, in 
computing Gf the largest error comes from Mw measurement with current value given as 
M w = 80.39 ± .06 GeV H and a much smaller error comes from Mz = 91.1867 ± .0020 



GeV. Although the error in W mass measurement is only 0.075% it gets amplified in 
Eq.© to 0.36% and thus constitutes the biggest uncertainty. In contrast, the error in 
Mz measurement gives an error of 0.015% in G F . From the above we can obtain the 
error in the determination of G F M to be 

G S F M = (1.166 ± 0.005) x 10~ 5 GeV~ 2 . (7) 

In the above equations the LFV physics and the SM are assumed to be factorizable. This 
is a good approximation since the S~ scalar does not couple to the W, because it is an 
SU(2) singlet, and hence does not affect the correction to the W propagator. We note 
that it does alter the running of a, but the effect is much smaller than the uncertainty 
in Mw and can be neglected. To get a bound on the couplings fu we demand that the 
corrections be no bigger the SM error given above. We use Eq. (P to obtain a bound on 
I/12I 2 as a function of Mj which is displayed in Fig.(jl]). 

Besides r M , the electron spectrum from muon decay also gives information on possible 
new physics. This is usually quoted in terms of the Michel parameter measurements. It is 
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interesting to note that only the following three Michel parameters get a correction from 
the Zee model 
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(9) 



(10) 



(11) 
(12) 
(13) 



All the others take their canonical SM values. The left hand sides of Eq.(|TT] - |T^) are 
experimental bounds taken from |TTJ . One other noteworthy point is that at the four-fermi 
level a right-handed current coupling is induced because of the S and vr interactions in 
lepton sector. However, this model does not have a corresponding coupling in the quark 
sector. This leads to a possible interesting scenario in which an apparent contradiction 
can occur in the positive outcome of right-handed current searches in the purely lepton 
sector such as in \i and r decays and a negative outcome in similar searches using hadrons 
such as in (3 decays of nuclei or pion decays. In Fig. [I] we also display the bounds on the 
couplings of the S boson from Michel parameters. These bounds are complementary to 
the limit one obtains from the lifetime study. 



Another tree level process that probes lepton number violation is r leptonic decay. 
The relevant effective four-fermi Lagrangian is obtained from Eq.(Q) by substituting the 
subscripts 2 to 3 and 1 to 3 appropriately. Here we use the branching ratio R e ^ defined 
by r — > euu/r — > \wv and obtain [compare Eq.(0)] 



M^(|/ 13 | 2 -|/ 23 
! 9 2 MI 

Since experimentally R eil is close to unity the error ±0.006 |L6[ only gives the absolute 
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value of the difference between I/13I 2 and | 1 2 - Notice that there is no dependence on 
the strength of the couplings involving u R . Thus, r leptonic decays can be used to probe 
a different region of parameter space from \i decays. Moreover, Eq.(|i~4|) shows that the 
/13 and are of the same order of magnitude. We also find the Michel parameter 
measurements here place less stringent bounds than from \x decays. 

Next we discuss the one loop effect involving the S-scalar. The most accurately mea- 
sured LFV neutral current process is the decay fi — > c-f. The branching ratio is calculated 
to be 



«(W| 2 +|/l3/23| 2 ) (My, , 

Br(/i ^ 67) = 2V KM 3 ) ' (15) 

This is sensitive to | <7i (72 1 2 + I/23/13I 2 and provides the most stringent limit [|DJ] on these 
parameters; see Fig. |2|. Since this is a sum of two positive terms and implies that \g\g2\ 
and I/13/23I must be of order 0.01 for M s = 800 GeV. We note in passing that similar r 
decays will test different combinations of coupling constants. 

A similar calculation gives the correction to the anomalous magnetic moment of the 
charged leptons, ag where £ = e,fi, or r. Then 

ae = ^ll£pEl(^±) 2 . (i 6 ) 



967T 2 \M S y 

Both a e and a M are very well measured quantities. However, since scalar interactions 
flip chirality and because of dimensional consideration we obtain two powers of lepton 



mass in Eg . (jig) . These powers of lepton mass, when combined with the loop suppression 
factors make these limits less constraining than some of the other processes considered 
here. Some of the other processes considered here, such as the \x decay parameters and 
the right and left handed lepton coupling to the Z discussed below, set tight limits on 
I/12I 2 , I/13I 2 , and I 1 2 - Therefore, we can use Eq (|T6| ) to put limits on g\ and g 2 . The 
limits from measurements of the electron and muon are shown in Fig. |3|. If the error can 
be reduced by a factor of 20 on a M the limit obtained is the lower curve in Fig. [| This is 
the target precision of the Brookhaven experiment E821 WE . 



We proceed to discuss the effect of the Zee model on precision measurement at the 
Z-pole. The main effect is on the left-right asymmetry of the charged leptons and lepton 
flavor changing neutral current decays. To one loop order the correction to the Zll vertex, 
ST^ is given by 

<5IV = 5^l7^ l + tfflkTAR > ( 17 ) 



where 



ffflL= 9^(§ l/fl|2 ) (^4- in6+m ) 



and 



2/^2 
W 



gfbs 



6 = _* £ ^_ (19) 
yK 288vr 2 v ; 

In the above we have used b = M^/Mg, sw = singly and cw = cos^iy- As expected 
the left-handed coupling is corrected in the Zee model without vr. However, if ur exists 
it will also modify the right-handed electron coupling to the Z in spite of the fact that, 
although u R is light, there is no tree level Zv c R v R coupling. Using Eqs. (p~8|) and (|T9|) one 
can calculate the left-right asymmetry, A LR . It is given by 

" Al H + (l-4^)[l + (l-4^) 2 ] 1 ' (20) 



where A^ R is the SM left-right asymmetry. The most stringent bounds are obtained by 



combined LEP [16|and SLC measurements [19| given in terms of the axial and vector 
couplings. This is shown in Fig. [|. 

With the above limits on the various couplings we can also predict the upper limits for 
the rare decays Z — ► lilj where U 7^ lj and U = e, fi, r. The width is calculated to be 

g 2 M z b 2 



24tt 5 (288c w ) 



\fikfjk\ (3\nb-c w ) + 



7T 



I I 2 4 

\9i9j\ s w 



(21) 



where i 7^ j 7^ k. With the limits given by the above considerations the upper limits for 
the e/z, fir, and er are all unmeasurably small in the foreseeable future. For example the 
upper limit for the Z — > efj, branching ratio is of order 10 -10 for a 800 GeV bilepton scalar. 

Having examined the constraints on the bilepton scalar we proceed to look into future 
tests. Short of discovering the S scalar itself, which is best done in a linear e + e~ collider, 
we propose that e, /1 ,r universality for the W boson will be a good place to look this kind 
of LFV physics. Explicitly, the branching ratios W — > evjW —>■ fiu and W —>■ eu/W — * tv 
will be very valuable. In the SM these branching ratios are unity at the tree level but 
at the level of first order radiative correction this universality is broken. However, this 
breaking is suppressed by the factor of |3D| a{m 2 /My^) or a(m%./Myy) and hence very 
small. The important point here is that the radiative corrections are accurately predicted 
in the SM, and the measurements of these branching ratios are very clean theoretical 
probes of LFV physics. We illustrate this in the Zee model by calculating the correction 
to the W leptonic decay widths. We find 



T(w -> lm) 



g 2 M } 



48vr 
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1 (— - II" 

2vr \ 3 12 



M 2 W 



72n 2 Ml 



1-31*4' 

Ml, 



(22) 



where i — e,fi ,or r. The first line is the radiatively corrected W boson leptonic width 



[21] and the second term in the curly bracket is the Zee model correction. As expected 



the left-handed nature of the coupling is preserved and v R does not play a role. We have 
also kept only the leading term in M^/Mj in the calculation and neglected the lepton 
masses. Eq. (^) immediately gives rise to the following results for the leptonic branching 
ratios: 

l + k(|f 23 | 2 -|f 13 | 2 ), (23) 
l + k(|f 23 | 2 -|f 12 | 2 ), (24) 

l + k(|f 13 | 2 -|fi 2 | 2 ), (25) 
where k = + 3 ln J^J ■ As 

seen from r universality the combination |/ 23 | — 
|/i 2 | 2 is bounded to be small and so we determine that the first branching ratio of Eq. (|23| ) 
can only accommodate a correction of less than 2xl0~ 5 for Ms= 800 GeV. Furthermore 
the constraints we obtained on |/i 3 | 2 and |/ 23 | 2 from the Z pole allow the universality 
violation in the last two branching ratios of Eq. (|2~3"D to be as large as lxl0~ 3 . At the 
LHC and the NLC where large samples of W decays are expected, these decays will be 
the most important and cleanest tool for probing LFV in the charged current sector. We 
note that the uncertainty due to the t-quark that plagues Ar in the interpretation of the 
muon lifetime measurement does not enter here at the one loop level. Once sufficient 
statistics are obtained these measurements will supersede many low energy tests of LFV. 

The last stop in our discussion of futuristic experiments is the production of the S- 
scalar. In particular we focus on pair production in e + e~ colliders. This is very similar 
to charged Higgs production in the two doublets model Since the Yukawa couplings 
fi2 and /13 are seen to be small the t-channel process involving neutrino exchanges can 
be neglected and one needs only to consider the s-channel virtual photon and Z exchange 
graphs. The coupling are all determined by the SM charges [see Table I] the production 
cross section can be calculated with Ms as the only free parameter. Explicitly, the cross 
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section for e + e — > S + S is 

na 2 f3 3 



a 



3s 



s(-l + 4s 2 w ) s 2 (-l+4s 2 w + 8s 4 w ) 
2c 2 w (s-M 2 ) + AcUs-M 2 ) 2 



(26) 



where s is the cm energy squared and ft — Jl — 4M|/s. If the mixing of the bilepton 
scalar with charged Higgs bosons can be neglected its dominant decay mode would be 
into liPj. The signal will be r and fi plus missing energy and unmistakable. 

As seen from the results we have presented, a bilepton scalar induces many interesting 
and testable effects. We have employed the Zee model augmented with a light i/r for 
quantitative studies since it is relatively simple and the number of free parameters are 
relatively few. We find that the strongest bound on the Yukawa couplings in the model 
comes from \x — > ej decay, lepton universality in r decays and the lifetime of the muon. In 
particular /12 and /13 are both constrained to be less than of order 0.1 and 1 respectively 
for Ms =800 GeV. However, the couplings gi have much looser bounds. Since Yukawa 
couplings are in general not universal we expect leptonic universality to be violated in W 
decays (see Eq. (0)). This illustrates the importance of W decays in probing LFV physics. 
Another effect which will be harder to discern is the modification of the U(l) coupling 
constant running. The S'-scalar contributes an amount of ^2 to the U(l) (3 function where 
g' is the U(l) gauge coupling. This will upset the unification of the SM gauge couplings 
at very high energies which is usually taken to be a hint for supersymmetry. We do 
not consider this to be a serious impediment for the following reason. Although the Zee 
model is interesting in its own right, it carries the same arbitrariness as the Higgs sector 
of the SM. Hence, we expect it be part of a larger structure that entails supersymmetry. 
Promoting the bilepton scalar to a superfield necessitates the introduction of a second 
scalar in order to cancel the anomaly caused by the fermionic partner of the S. In short 
one would have to enlarge the minimal supersymmetric standard model by two bilepton 
superfields and also one singlet neutrino superfield. The details of how this can be achieved 
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and the ensuing intricate phenomenology is beyond the scope of the present paper and 
we shall defer the study of this issue to a later work. 

This work has been supported by the Natural Sciences and Engineering Research Coun- 
cil of Canada. 
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Table 1: Lepton number, SU(2), and U(l) charges of the left-handed lepton doublet, the 
vr, the charged lepton conjugate, and the S-scalar 
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Figure 1: Bounds on coupling constants derived from measurements of Michel parameters 
for muon decay are plotted as a function of the mass of the scalar, Ms- The constraint 
on I/12 1 2 comes from the measurement of the muon lifetime. 



Figure 2: As in Figure 1, for the constants I/23/13I 2 + \giQ2\ 2 - This constraint comes from 
the experimental limit on the width of the decay process \i — > cy. 



Figure 3: Coupling constants plotted as in Figure 1. These constraints come from precision 
measurements of the anomalous magnetic moment of the electron (top line), and of the 
muon (middle line). The lower line shows the bound which could be obtained if the muon 
measurement were twenty times more precise. 



Figure 4: Coupling constants plotted as in Figure 1. This plot represents three different 
constraints as determined by the choice of leptons, e,/z, r for the parameters, l,j, k (with 
I 7^ j 7^ k). These constraints are derived from the measurements of leptonic axial and 
vector couplings. 
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